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►^ . Abstract 

1^ \ Coexisting periodic solutions of a dynamical system describing nonlinear 

O ■ optical processes of the second-order are studied. The analytical results con- 

^r , cern both the simplified autonomous model and the extended nonautonomous 

vQ \ model, including the pump and damping mechanism. The nonlinearity in the 

r^ ' coexisting solutions of the autonomous system is in concealed frequencies de- 

pending on the initial conditions. In the solutions of the nonautonomous 
system the nonlinearity is convoluted in amplitudes. The neighbourhood of 
^ ' periodic solutions is studied numerically, mainly in phase portraits. As a re- 

K^ ■ suit of disturbance, for example detuning, the periodic solutions are shown to 

K> , escape to other states, periodic, quasiperiodic (beats) or chaotic. The chaotic 

^ \ behavior is indicated by the Lypunov exponents. We also investigate selected 

aspects of synchronization (unidirectional or mutual) of two identical systems 
being in two different coexisting states. The effects of quenching of oscilla- 
tions are shown. In the autonomous system the quenching is caused by a 
change in frequency, whereas in the nonautonomous one by a change in am- 
plitude. The quenching seems very promising for design of some advanced 
signal processing. 
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1 Introduction 

Nonlinear systems are usually characterized by two or more coexisting states, corre- 
sponding to the same values of parameters. As for the first time found by Poincare', 
a periodic solution disappears (or appears) by couples in the form of real roots of an 
algebraic equation. Today it is well known that both chaotic and periodic states can 
exist within standard nonlinear models like Duffing oscillator P , Lorenz and Rossler 
systems or Henon map j2]. The dynamical coexistence is frequently referred to as a 
generalized multistability [HI E] . Multistable behavior appears also in nonlinear op- 
tics jS], electronic circuits mechanical systems [7j and neurons nets [Hlini- There 
is no universal (global) method for detecting multistability of a dynamical system. 
Usually, we do this numerically, for example, looking for different basins of attraction 
in a system with selected parameters. Obviously, the numerical approach is neces- 
sary if we try to find coexisting chaotic states. In some cases, however, analytical 
form of coexisting states (solutions) is also possible to find, provided that the states 
are regular (periodic). Frequently, purely numerical investigation does not deliver 
sufficient information about the intricate nature of coexistence. Therefore, attempts 
at finding some analytical results, if it is possible, are always physically valuable. 
We try to follow such an approach in this paper. The aim of this study is to find 
a class of coexisting periodic solutions in a well known nonlinear model describing 
the nonlinear optical processes of the second order and to investigate the behaviour 
of these analytic solutions in response to a disturbance of the source of periodicity. 
The dynamical model is considered in a simplified version, that is autonomous and 
nonautunomous ones. The dynamics of the autonomous and nonautonomous sys- 
tems are compared within the phase space. The coexisting periodic solutions in the 
nonautonomous model are shown to be controlled by changing the system's param- 
eters, which leads to transitions from one state to another. The coexisting states 
can escape to chaotic or nonchaotic (periodic) states. The type of the final state 
is deducted from the Lypunov exponents. Finally, the behaviour of two identical 
dynamical systems being in two different coexisting states on a linear interaction 
between the systems turned on, is studied. In some cases one or two of the coexist- 
ing states can be quenched. The quenching effects are shown to be controllable by 
the parameters of the system. 

2 Equations of motion 

Let us consider a nonautonomous dynamical system governed by the following set 
of equations [IDl EH Ca CH! : 

— = -iooa - -fia + ea*b + Fie~'^^\ (1) 

^ = -^2ub - ^2b - ^ea' + F^e^'""-'. (2) 



Physically, the equations describe an interaction between two optical modes of the 
frequencies uj and 2uj. The complex dynamical variables a and h are the ampli- 
tudes of the fundamental and second-harmonics modes, respectively. The interac- 
tion takes place via a nonlinear crystal placed within a Fabry-Pefot interferometr. 
The quantity e is a nonlinear coupling coefficient, whose value is proportional to the 
second-order nonlinear susceptibility. The parameters 71 and 72 are the damping 
constants of the fundamental and second-harmonics modes, respectively. Moreover, 
the system is pumped by two external fields Fic"**^^* and ^26"*^^*, where Fi and F2 
are electric field amplitudes at the frequencies Vti and 0.2- respectively. Henceforth, 
all the parameters, that is tu, e, Fi^2, and VLi^2 are taken to be real as in |TTj. 
To visualize the dynamics of the system(^-(j21) the four-dimensional space 
{Re a. Re b, Im a, Im b) 

is required; but, as it is impossibile, we carry out the visualisation for its two di- 
mensional sections. The system(PQ)-(j2I) does not belong to the class of integrable 
systems and usually it is studied numerically. However, in special cases, analytical 
solutions of (UJ-® are also possible. Below, we consider a class of coexisting peri- 
odic solutions of the system and qualify the kind of motion in their neighbourhood. 

2.1 Autonomous case 

Let us first consider the problem of periodic orbits in the simplest (conservative) 
version of the system(|H)-(P)). that is 

— = —tuia + ea b , [6) 

db , 1 9 

- = ^r2.b--.aK (4) 

The equations of motion Q^Q were used for the first time by Bloembergen to 
describe second-harmonic generation of light ^2 El E] • The above system has two 
coexisting periodic solutions (the details may be found in Appendix A). The first 

a«(t) = ae-*("+^^^'^^)* , (5) 

5(i)(t) = _'_J?.e-^^i^+h^/^)^ (6) 

and the second 

a(2)(t) = ae-*(— l^v^^)*, (7) 

6(2)(t) = +lJ^^e~'^^^-¥^/^)t_ (8) 

The coexisting first harmonics © and ((Zj) have identical amplitudes a and different 
frequencies u + 0.5ev^a*a and oj — 0.5e-\/a*a, being functions of the initial condition 
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Figure 1: Phase-diagrams of the coexisting solutions ©-(0) and ((Tj)- (jHl) in black 
and red, respectively. The parameters are: i?e a = 10, /ma = 0, cj = 1 and e = 0.1. 
The positions of the phase points 1 and 2 are marked at the initial time t = and 
at the later times t = 0.8 or t = 2.5. 
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The second harmonics © and (jH)) have the same amplitudes 



and different 

coexisting frequencies 2a; + e\/a*a and 2uj — e\/a*a. Additionally, the function 6^^^ 
is of the sign opposite to that of b^^^ . 

Let us note that the functions ((Tj) - (jH)) can be constant in time (the period T = oo) 
ii Lj = ^e\/a*a. Physically, it means that the vibrations are quenched. The fact 
that a frequency (period) depends on amplitude is well known in the theory of au- 
tonomous systems J7j . A variation of the period with amplitude is well known, for 
example, in the case of a pendulum for larger deviations. 

In the four- dimensional phase space {Re a, Re b, Im a, Im b) the solutions (0)-® 
and ©-(jHl) generate two coexisting hyper- surfaces. The geometrical relationship 
between them are illustrated in six two-dimensional phase diagrams (Fig^, where 
the coexisting solutions create simple Lissajous-like curves. Some of them are iden- 
tical (degenerate), as readily seen, for example, in Fig. la. Both phase points 1 and 
2 start together from the same position (t = 0), rotate in the same direction and 
draw the identical orbits {Re a^^^)^ + {Im a^^^'Y = r^ and {Re a^'^^Y + (-^"^ a^'^^Y = r^, 
where r^ = a*a. The only difference is, that the phase point 1 draws the circle at 
the frequency u + ^ey/a*a, whereas the phase point 2 at the frequencycj — ^e^/a*a. 
The other Lissajous-like trajectories are presented in Fig. lb - Fig. If. As seen in 
Figs.(a,c,e), points 1 and 2 start from the same position, whereas in Figs.( b,d,f) 
from two different ones. 
Generally, if we start from any point lying on a periodic trajectory of the sys- 



tern Q- Q we always remain in the same trajectory. This is frequently called the 
translation properties of autonomous systems which means that to a given trajec- 
tory corresponds an infinity of motions (solutions) differing from each other by the 
phase p/^. The question is, however the type of trajectories of the system dS} - 
(^ when the initial conditions do not lie on periodic orbits. There is no difficulty 
in solving this problem numerically in contradistinction to a general analytical ap- 
proach. The numerical analysis shows that all trajectories of the system (jH)) - (P)) 
originating from independent initial conditions a(0) = a and 6(0) = /3 (contrary to 

the dependent conditions a'^^'^)(0) = a and 6(^'^)(0) = ±W^ generating periodic 
orbits) behave in a nonperiodic, mainly quasiperiodc manner. For example, when 
a(0) = a and 6(0) = we get an analytical result known since the pioneering work 
by Bloembergen ^1 Ej : 




a{t) = a seek -^aet e""^* (9) 



h{t) = — ^ tanh l—aetj e"'""* (10) 

where a = a*. Obviously, due to the time dependent amplitudes asechi^aetj 

and — -% tanh f ^aet) the functions a{t) and b(t) are not formally periodic i.e. 
a{t + T) 7^ a{t) and b{t + 2T) ^ b{t), where T = co'/27r. The solutions are frequently 
referred to as nearly periodic, which mearly means that in the course of time both 
functions approach a purely periodic motion. 

Another question is the change in the harmonic solutions (jS))^® on inclusion of 
damping in the dynamical system ©^Q: 

da . ^, ,^^, 

— - = —tuja — 'ya + eab, (11) 

dt 

db , , 1 9 , , 

_ = ^^2u;b - jb - -ea\ (12) 

where, for the sake of simplicity, we have put 71 = 72 = 7- Using the method 
[T7j proposed by Krylov-Bogolubov (we do not present the details here) we get two 
coexisting nonperiodic solutions {a^^\t) , b^^'' (t)} and {a^'^\t) , b^'^\t)} , namely: 

^ ^ 2\ a* ^ ^ 

In the limit 7^0, the above functions become periodic solution (0)-®- Here, in 
contradistinction to the periodic solutions, not only the amplitudes are functions of 
the damping constant 7 but also the phases. For t —>■ 00 the coexisting solutions 
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Figure 2: Phase diagram of the function ()T3jl for 7 = (circle) and 7 = 0.2 (spirals). 
The circle is identical to that in Fig. la. The phase points 1 and 2 escape from the 
joint coexisting orbit due to damping. The first (black) goes faster then the second 
(red). 



{a(^)(t),6(^)(t)} as well as {a(2)(t), 6(2)(t)} tend to zero. 

Damping damages degeneration of the periodic orbit in Fign(a,b,c and e). By way 

of example, it is illustrated in Fig. 2. For 7 7^ the phase points 1 and 2, drawing 

different trajectories (1 goes faster than 2), approach the fixed point (0, 0), being an 

attractor. 

Equations (fTT|) - (fT^ may be solved subject to the initial conditions a(0) = a and 

6(0) = to yield 



a{t) 



hit) 



ae '''* sech I - — aeil 



ae 



-It 



^ 




-■yt\ 



.,—iujt 



tanh I ^^aefl - e-T*) e"*^-'* 



(15) 
(16) 



The behavior of the above amplitudes and phases is remarkably different from that 
presented by Eqs. (fT^ - (fT^ . Here, the phases depend neither on the initial conditions 
nor on the damping constant. The amplitudes are damped in a much more intricate 
way than those in Eqs. (fT ^ -(fT ^ . 

In the phase space, both functions (fT^ -(fT ^ and (P3|) -(fTH |) tend to the same attractor 
being a fixed point. Obviously, in the limit 7 — >■ the functions ()15|) - (fT^ approach 
arbitrarily closely functions (Pl- ffTUI) . respectively. Generally, the numerical studies 
show clearly that in the phase space the system (fTT|) - (fT^ always tends to a fixed 
point, independently of the initial conditions. 



2.2 Periodic resonance solutions of nonautonomous system 

In order to find a periodic resonance solution of the systeni(^- (J2|) we look for a 
solution in the form 



a{t) = Xe 



-iQ-\t 



b{t) = Ye 



-iD,2t 



(17) 



where VLi= uj and 0.2 = 2u are the resonance conditions. The amplitudes X and Y 
are constant in time. On inserting (|T7jl into (jT))- Q we get two algebraic equations 
(quadratic) in the complex variables 



- 7iX + eX*Y + Fi = , 



-72^ - -tX' 



0. 



(18) 
(19) 



Therefore, we look for X and Y as functions of the parameters: e, 71, 72, -Fi and F2. 
Restricting the number of the parameter to three, we get solutions whose physical 
context is clear, and the algebraic form is easy for numerical investigation. 



2.2.1 The case I, Fi = , F2 ^ , 71 ^ and 72 ^ 

This case {Fi = 0) describes the subharmonic generation. There are two coexisting 
solutions {a'-^\t),b{t)} and {a^^\t) , b{t)} 



a(i'2)(t) 



T\ 



I2F2 27172 



-iujt 



Fo> 



7i72 



h{t) = ^e 
e 



7l ^-i2ujt 



(20) 
(21) 



which differ only in the phase a^'^\t) = a^^\t)e ™ . Therefore, the coexistence has a 
trivial character, (see Eqs.(3.1, second line) in Ref.|lUj). 



2.2.2 The case II, Fi 7^ , F2 = 0, 71 ^ and 72 ^ 

Physically, this case {F2 = 0) corresponds to the second harmonic generation. There 
are no coexisting solutions but only one single solution (see Eqs.(3) and (4) in 
Ref.HH): 



a{t) 
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V e2 



(22) 
(23) 



2.2.3 The case III, Fi ^ , F2 ^ 0, 71 ^ and 72 = 

Here, both the subharmonic effect and the second harmonic processes compete with 
each other. The assumption 72 = describes the so-called good frequency conversion 
limit for subharmonic generation. It is easy to prove that the system (Q)-(j21) has 
two coexisting solutions {a^^\t) , b^^\t)} and {a^'^^t) , b^'^\t)} given by 

«"*(() = T^e-"', (24) 

Physically, for the same values of parameters the system has two periodic states 
differing in the vales of amplitude. If ^ = ^ y ^; then second-harmonic vibrations 
are quenched (b^'^^t) = 0). In this case the subharmonic generation is maximal. 

Phase diagrams for the coexisting solutions {a^^\t) , b^^\t)} and {a^'^\t) , b^'^\t)} 
are presented in Fig. 0. As seen, only the phase curve for the pair {Rea^^\t), Ima^^'{t)) 
covers the phase curve for the pair {Re a^'^\t) , Im a^'^\t)) . Here, point 1 is by vr out 
of phase with point 2. The other curves are non-degenerate that is they are sepa- 
rable. This behaviour follows from the fact that the functions b^^''^\t) differ in the 
values of amplitudes. 
There is no geometrical correspondence between the phase portraits in Fig^ and 

FigH 

The main difference is, however, between the autonomous and nonautonomous phase 
dynamics. The translation properties of autonomous systems ( sometimes called free 
phase PHI) do not hold in nonautonomous systems. It means, for example that the 
phase points 1 and 2 in FigI21 (a) follow the circle provided that they start only from 
the points (a(0) = -10, 6(0) = 0) and (a(0) = +10, 6(0) = 0), respectively. If they 
start from the other points lying on the circle they escape from it, which does not 
take place in the autonomous case (this problem is considered in detail in Section 
HI). 

As seen from Eqs. (j24j) -(|25| ) the parameter e governing the nonlinearity of the sys- 
tem (Q)-© is felt in the amplitudes, in contradistinction to the autonomous case, 
where e is felt in the phases (see Eqs.flT ^ - (fT^ ). 

2.3 Fractional resonance 

It is assumed that the difference between the periodic solutions of the autonomous 
and the nonautonomous system is mainly that the solution of the former has the 
period (frequency) being a function of the initial conditions and the parameter gov- 
erning the nonlinearity of the system itself (vide Eqs.©-®), whereas that of the 
latter have period of the external pump fields only (vide Eqs. (l24j) - (f^H|) ).The main 
difference between the periodic solutions of the autonomous and nonautonomous 
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Figure 3: Phase-diagrams of the coexisting solutions (f^ - (f^ for cu = 1, e = 0.1, 
7i = 0.002, Fi = 0.01 and Fa = 5. 

systems is that the period of the former is a function of initial conditions and the 
parameter governing the nonlinearity of the system (vide Eqs.dnj)-®), while that 
of the latter is determined by the external pump fields only (vide Eqs. ((211) -(EH))- 
However, in some cases the periodic solution of the nonautonomous system may have 
the period dependent also on the parameter governing the nonlinearity of the dif- 
ferential equation. This takes place in a special case of the resonance, namely when 
we want a nonautonomous system to vibrate at the frequency of its autonomous 
counterpart. Then, instead of looking for the solutions of (^Q)^© in the form of (fT7|) 
we search for solutions given by (see Eqs.®-®) 



a{t) = xe~''''\ b{t) = T-. 



-i2nit 



(26) 



where Qi = u ± \e^x*x. It is easy to note that the functions (|^ satisfy Eqs.(IH)- 
provided that 72 = 0, Fj = and x = Fi/71. In this way the frequency fli 
becomes additionally a function of the damping constant 71 and the amplitude Fi. 
By way of example, for 72 = 0, F2 = 0, cu = 1, e = 0.1 and Fi = 5 the system ((H)- 
^ has a periodic solution provided that fii = 1 ± |. Therefore, we have two sets 
equations. The first 



-^ = -ia - 0.5a + 0.1a*b + 5e-'^\ 
at 



db 
dt 



-i2b - 0.5ea^ 



(27) 
(28) 



where the solutions are given by a{t) = 10exp(— i|t) and b(t) = — 5iexp(— Sit) and, 
the second 

^ = -ia - 0.5a + 0.1a* b + 5e-''^\ (29) 

— = _^26-0.5ea^ (30) 

(Jjb 

where a{t) = 10exp(— i|t) and b{t) = 5iexp{—it). Both sets of equations describe 
the so-called fractional (subharmonics, demultiplication) resonances. The periodic 
solutions of (J77j) - (PH|) and ((221)^(1201) satisfy the conservative autonomous system 
da/dt = —ia + 0.1a*b and db/dt = —i2b — 0.5a^ and have phase representation 
identical to that in Fig.((T)). Finally let us note that the system (III)-(|2I) if 72 = 0, 
F2 = and the resonance condition fii = uj holds, has a Bloembergen-type solution 
in the form: 

a{t) = ^sech(——et]e-'^\ (31) 

71 V 2 71 / 

Kt) = -^tanhf^^et)e---. (32) 

7iV2 V 2 71 y 

The above functions for t —>■ 00 tend to periodic states. 

3 Chaotic behaviour 

The coexisting periodic solutions P^ - (j^ naturally lead to the question about the 
effects of a disturbance of the resonance conditions. It is intuitively clear that this 
problem can only be solved numerically. As a background to numerical investigations 
we use the equations 

-^ = -m- 0.002a + 0.1a*6 + 0.01e"^^^*, (33) 

at 

^ = -226 - O.OSa^ + 5e-*^2* , (34) 

dt 

where Qi and ^2 play a role of parameters. If at the time t = the state of the 
above system is determined by the initial conditions a(0) = 10 and 6(0) = 0.01 and 
the conditions of resonance are satisfied i.e. fii = 1 and ^^2 = 2, then the pair of 
periodic functions 

a{t) = lOe"^* , 6(t) = O.Ole-^2*. (35) 

satisfy the differential equations ((221)^ (EH)- Phase diagrams of the functions (jH^j) are 
given in Fig.3 (black lines). For the initial conditions a(0) = —10 and 6(0) = 0.03 
we get the second pair of coexisting periodic solutions (red line in Fig.3) 

a{t) = -lOe"^* , 6(t) = 0.03e-^2* (36) 

10 



Let us now consider the behaviour of the system in the neighbourhood of the reso- 
nance periodic solutions in the range 

fii = 1 and < ^2 < 4 . (37) 

Numerically, we solve Eqs. p3|) - (J3l|l . with the help of a fourth-order Runge-Kutta 
method, with the initial conditions (a(0) = 10, 6(0) = 0.01), for fii = 1 and selected 
values of ^2 7^ 2. Then, we repeat the same procedure with the initial conditions 
(a(0) = 10, 6(0) = 0.03). The numerical results are reflected in phase portraits and 
in the spectra of Lyapunov exponents, computed by the Wolf procedure jTH] • 

At the beginning, let us consider two the most typical types of behaviour being 
a result of detuning, that is the resonance condition breaking. 

Example I - large frequency detuning. 

Take the frequency ^2 = 0.472 instead of Q2 = 2. We observe that initially the phase 
point follows the periodic solution (jH^ and then escapes from it to a chaotic state. 
This result of the large detuning is seen in the phase {Im a , Rea) plane FigEJ^a). A 
symmetric behaviour is observed in the same phase portrait, if the system (J33|l - P^ 
starts with the initial conditions a(0) = —10 and 6(0) = 0.03, see Fig^b). The peri- 
odic (^2 = 2) and chaotic {Q2 = 0.472) states of the system presented in FigE] con- 
firm the spectra of the Lyapunov exponents{-0.0007, -0.0007, -0.0022, -0.0022} 
and {0.1659,0.0012, —0.0029, —0.1699}, respectively. The latter spectrum, contain- 
ing two positive exponents, indicates a strong chaotic behaviour, the so-called hy- 
perchaos. 

Example II - weak frequency detuning 

Take now the frequency Q2 = 2.25 instead of Q2 = 2. Here, the point escapes from 
the periodic state (governed by Eqs. ljH^jl Ho another periodic state characterized by 
the spectrum {-0.0013, -0.0014, -0.0015, -0.0015}. This is shown in (FigEj). The 
same phase structure is obtained if the system starts from the coexisting periodic 
states described by Eqs. (jnU|) ). 

Globally, the behaviour of the system in the resonance neighbourhood ()37|) is pre- 
sented in Figini The spectra of the Lyapunov exponents {Ai, A2, A3, A4} versus ^2 
show the regions of order or chaos. If Ai > then system is chaotic (black colour), 
if simultaneously Ai > and A2 > then system is hyperchaotic (black and green) 
and finally if Ai < the system behaves nonchaotically (periodically). For ^2 = 2 
the system is in a periodic state fEqs. ljH^ or solid lines in Figl^I). 
Fig. (jHl) shows that the spectra are nearly symmetric relative to the resonance fre- 
quency Q2 = 2. By way of example, for ^2 = 2.25 (see Figl^I) we have the same 
spectrum as for ^2 = 1.75 and get two identical phase structures. 

A convenient way of finding out what may be expected in the system ()33|) - (jMj) 
under conditions fii 7^ 1 and f22 7^ 2 is to calculate the maximal Lyapunow exponent 
Ai as a function of Qi and ^2- Then, the global dynamics of the system is simply 
presented by the Lyapunov map in the space of (r2i,r22) (Fig|7j), where the values 
of Ai are marked by an appropriate colour. A band structure of the map shows that 
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Figure 4: The symmetric escape of the phase points from the coexisting orbits. 
Figure (a) - evolution of the system (|H ^ - (jn^ with the initial conditions a(0) = 10, 
6(0) = 0.01 and for 1) i7i = 1, f22 = 2 (periodic solution - read) and 2) Qi = 1, 
Q2 = 0.472 (chaos - black). Figure (b)- the same as in Figure (a) but with the initial 
conditions a(0) = -10 and 6(0) = 0.03. 




Figure 5: Escape from one periodic state (red) to another periodic state (black). 
Evolution of the system (IHH jl - ljH^ with the initial conditions a(0) = 10 and 6(0) = 
0.01 and for 1) fii = 1, ^2 = 2 (red) and 2) ^i = 1, ^2 = 2.25 (black). Black orbit 
has been drawn for t > 4500, to omit the transient effects). 
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Figure 6: (a) Spectrum of Lyapunov exponents Aj, i = 1,2,3,4, for the system 
dMI-dSll) with the initial conditions a(0) = 10 and 6(0) = 0.01, for fii = 1, and 
< 1^2 < 4. (b) An enlargement of the region of | Aj |< 5 x 10^'^. 



the system is much more sensitive to the changes in Q2 than to those in Qi. This 
is caused by the fact that the amplitude of the forcing term in (jH^ is 500 times 
lower than in the forcing term in (jS^)- Therefore, the changes in Qi do not affect 
so much the stability of the system as those in 0,2- The map is simply symmetric 
to the line f^i = 1 which means that the detuning Q2 + S as well as ^^2 — 5 causes 
the same kind of stability (instability) of the system. The cental point of the map 
{Qi = 1,1)2 = 2) corresponds to the periodic solution (jH^j) . 

We might expect that if the nonlinear interaction is sufficiently weak, that is when 
e/u << 1, then in the presence of detuning (that is in the neighbourhood of the 
periodic solution (jHS))) beats could appear. Though the nonlinear interaction in (j^ - 
(jM|) can really be treated as weak (u; = 1, e = 0.1) it is not sufficiently weak to induce 
quasiperiodic phenomenon. Distinct beats appear in (|H ^ - (jn^ when e/u < 0.01. 
This condition is satisfied (with a wide margin) if we put in ()33 |) -()34 |) . for example: 
Qi = uj = 100 (instead oi u = 1) and Q2 = 200 ± 6, where 6 is the value of detuning 
(for 6 = the system has the periodic solutions redefined a{t) = 10e~*^°°* and 
b(t) = 0.01e~*^°°*). A typical example of beats is presented in FiglHl The beats 
proceed as follows: if ^2 tends to 2uj = 200 i.e. {6 -^ 0) the amplitude of beats 
increases but this happens only to a certain value of the difference 12^7 — 021 = A 
(in our case A = 0.5) after which the amplitude of beats began to decrease. Finally 
beats disappear for 2uj = Q2 = 200 that is in the resonance. The beats problem 
(also chaotic beats) in different nonlinear systems has been recently investigated in 
nonlinear optics [201 ^^^ i^ electric circuits [21]. 

Nonautonomous systems do not manifest the so-called translation properties as 
their autonomous counterparts do. This property can be readily refiected in phase 
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Figure 7: The values of maximal Lyapunov exponents for Eqs. ()33|) - (jn^ with the 
conditions a(0) = 10 and 6(0) = 0.01 if < fi < 2 and < fia < 4. 
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Figure 8: Beats. Time evolution Reh{t) of the system (I!))-© where oj = 100, 
7i = 0.002, 72 = 0, Fi = 0.01, F2 = 5, fli = 100 and the initial conditions are 
a(0) = 10 and 6(0) = 0.01 and for (a) ^2 = 202, (b) (^2 = 200.1. The red band in 
(b) represents the periodic solution ^2 = 200. 
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Figure 9: Escape of the phase curve starting from the point B lying on the periodic 
curve (red circle). The periodic orbit, generated the system (jHH|) -(|H ^ . starts from 
the point A (the initial conditions: a(0) = 10 and 6(0) = 0.01). The same system 
on starting from B (the initial condition) a(0) = —0.0922 — ilO.0196 and 6(0) = 
0.0102 - «0.002 behaves chaotically. 

portraits. By way of example, to demonstrate this behaviour we use Fig. 3a (black). 
If the system §^- §^ starts from the point A = (a(0) = 10, 6(0) = 0.01) then in 
the phase portrait we observe a circle (the periodic solution ()35|)). However, if the 
same system starts from another point B ^ A, lying on the circle, it does not remain 
on the circle but escapes from it and draws another curve - in our case of chaotic 
one (FigEl). This chaotic behaviour is confirmed by the spectrum of the Lyapunov 
exponents {0.0059, -0.0009, -0.0031, -0.0077}. 

Generally, the system can escape from the periodic orbit to stable (periodic) states 
or to chaotic states. On calculating the maximal Lyapunov exponents Ai for the 
system starting from individual points of the periodic orbit we get the information 
on which orbit is chaotic Ai > and which is nonchaotic Ai < 0. This is shown in 
Fig.dinD. 

4 Synchronization of the coexisting states. 
Quenching 

Let us now consider the synchronization problem (mutual or unidirectional) of two 
dynamical systems (a, 6) and {A,B), which we may assume to be identical in all 
respects but being in two coexisting periodic states. We take the system (a, 6) given 
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Figure 10: Maximal Lyapunov exponents Ai > for the trajectories of system (jH!^ - 
(jM|) started from different points of the periodic orbit (red) in Fig. 9. The individual 
values of Ai show from which points (lying on the orbit) chaotic (Ai > 0) , and from 
which nonchaotic (Ai = 0) phase curves originate. 



by 



and its copy {A, B) and couple them linearly: 
da 



dt 

dh 
'dt 

dA 
It 

dB 

m 



-iLoa 



S: 



{a,A) ( 



7ia + ea*b + Fie"*^^* 



= -i2ujh - 726 - -ea^ + Fae'^^^* 

- 5(5,B)(6-5) 

= -iuA - 7iA + tA*B + Fie-*^^* 

- S(^A,a){.A- a) 

= -i2uB - 72^ - ]^eA^ + Fae"'^^* 



{3t 



(39) 



(40) 



(41) 



where ^(jj) = Sq-^,) is a parameter of the coupling. The couphng is usually turned 
on at an arbitrarily chosen time. 

The most spectacular behaviour is observed if the systems (a, b) and {A, B) are au- 
tonomous and conservative, that is when 71 = 72 = and Fi = F2 = 0. Suppose 
that the system (a, b) is in the state (0) - © whereas {A, B) is in the state ((7j) - ©• 
This means that the oscillator a vibrates at the frequency uj+5, where 5 = 0.5ey/a*a, 
whereas the oscillator A vibrates at the frequency 00 — 6. The oscillators b and B 
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Figure 11: The temporal oscillation death in a, A-components, black and red, re- 
spectively (a); and oscillation in 6, 5-coniponents (b). The effect forced by mu- 
tual synchronization in the system (|38p - (j4H) with the initial conditions a(0) = 10, 
6(0) = 5, A{0) = 10, 5(0) = -5, for cj = 1, e = 0.1, 71 = 72 = 0, Fi = Fa = and 
SaA = SbB = SAa = SBb = 0.5. The ^-terms are turned at t = 10. 



vibrate at the frequencies 2(u; + 6), 2{uj — 6), respectively. Therefore, the pairs(a, b) 
and {A, B) are simply detuned in frequencies. If the oscillator pairs (a, h) and (A, B) 
are mutually coupled, they get synchronized. This is illustrated in FiglTTl As seen, 
the coupling S = 0.5 is turned on at t = 10 and the systems synchronize at t = 15. 
Consequently, for t > 15 a(t) = A(t) and b(t) = B(t). What is more, the oscillator 
a and A are quenched for some time but at the same time the oscillators b and B 
vibrate periodically with the double frequency 2uj. As seen from Fig^J after the 
oscillation death in a, A-components we always observe its rapid and short revival 
- these effects occur one after the other. The revivals always correspond to appro- 
priate collapses in b, 5-components The quenching interval depends on the coupling 
constant S, the lager the value of S the longer the quenching interval. Therefore, the 
quenching effect is the most effective in the case of strong interaction that is when 
S » uj. The quenching phenomenon is forced by the real parts of the S'jj-terms in 
ISHIHTI They are simply the momentum (velocity) terms, sometimes named diffusive 
j22j, being sources of dissipation in the individual systems. 



Complete quenching can be obtained in the case of unidirectional interaction. 
Suppose, that S(B,b) = S(A,a) = 0, which means that {A, B) is a transmitter (master) 
system and (a, b) is a receiver (slave) system. Moreover, the transmitter sends a 
continuous signal [22] that it is in the state described by ((Tj) - (jHI) because we 
choose the frequency uj = ^e\/a*a. It is possible, for example, ii cu = 0.5, a = 10, 
Sa,A = Sh^B = 0.5 and e = 0.1 fFiglT^. Physically, it means that the transmitter 
does not vibrate. The receiver is in a periodic state © - (|S1)- On turning on the 
coupling the vibrations in the receiver are quenched - oscillation death is complete. 
The total quenching can be spectacularly observed with the help of the phase plots 
in Fign(a). Namely, the phase point 2 does not move, whereas point 1 follows the 
periodic orbit. At the moment we turn on the synchronization mechanism and point 
1 moves slower and slower towards point 2, and finally approaches the point 2 - the 
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Figure 12: The complete quenching caused by unidirectional synchronization in 
the system(jnHI)-(|lII) with the initial conditions a(0) = 10, 6(0) = 5, A{0) = 10, 



B{0) = -5, for io = 0.5, ( 
The terms S(a,A) = S(b,B) 



= 0.1, 71 = 72 = 0, F, = F2-- 
0.5 are turned on at t = 10. 
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Figure 13: Return to the periodic orbit of a point chaotically escaping from it. The 
escape is shown in FigEl The return is caused by unidirectional synchronization 
switched at the time t = 55. 



orbit is quenched. A similar behaviour is observed in FigH(b), (c) and (e). In FigQ 
(d) and (f ) point 2 does not lie on the orbit of point 1. Point 1 moves up and down 
along the parabole. If we turn on the synchronization mechanism point 1 escapes 
from its orbit and tends to point 2 - the orbit is quenched. Also, the inverse process 
is possible, that is creation of an orbit for point 2. 

With the help of the unidirectional synchronization we can force the phase point 
in FiglHlto return and escape chaotically from the periodic orbit. To do that, it is 
necessary to use the svstem (J38|) - p9|) with S(A,a) = S[B,b) = 0, 71 = 0.002 ,72 = 0, 
Fi = 0.01, F2 = 5, fii = 1, ^2 = 2 and the'initial conditions a(0) = -0.0922 - 
ilO.0196, 6(0) = 0.0102 - ^0.002 (chaotic orbit) and a(0) = 10 and 6(0) = 0.01 (red 
circle). Moreover, the terms S(a,A) = S(^b,B) = 0.5 should be turned on at t = 10, 
then we observe the situation presented in FiglT^ 

Synchronization of the nonautonomous coexisting states is as easy as that of the 
autonomous ones. The resultant vibrations usually appear as intricate revivals and 
collapses, frequently chaotic. 



5 Conclusions 

The coexisting periodic solutions, presented and analyzed in this paper, are a natural 
feature of the nonlinear differential equations (0)-(j21) describing nonlinear optical 
processes of the second order. 

The coexisting solutions of the autonomous type (0)-® have a simple physical 
interpretation. Namely, we are able to prepare the initial states of two independent 
harmonic oscillators (of the frequencies uj and 2uj) in such a way that after turning 
on the nonlinear interaction the oscillators vibrate at the detuned frequencies u + 5 
and 2{ijj + 5) or at the coexisting frequencies oj — 5 and 2{ijj — 5). If the system 
is damped, this correction also depends on the damping constant. This is clearly 
seen from Eqs. ()13p - (fT^ . The coexistence in frequency but not in amplitudes seems 
to be characteristic of a large class of nonlinear autonomous systems in the so- 
called rotating wave approximation (see Appendix). Switching on a linear coupling 
between the two coexisting detuned systems leads to the temporal disappearance 
of {uj ± 5)-vibrations and the appearance of 2c(j-vibrations only. The quenching 
interval can be controlled by the coupling parameter. In special cases (unidirectional 
coupling) we are able to completely annihilate vibrations in both oscillators being 
in the coexisting states. 

The physical interpretation of the coexisting solutions for the nonautonomous 
case is different to that presented above. The nonlinearity is now concealed in am- 
plitude, and the frequency is a parameter of the system. Two periodic solutions, 
having the same resonance frequency, have different amplitudes (or the same am- 
plitudes but different signs). Therefore, the coexistence means here a possibility of 
existence of two resonance states at the same values of parameters. The structure 
of the coexisting states in the phase space is completely degenerate as in the case of 
Eqs. fl20|) - (PT|) . The difference is only in phase, or partially degenerate as in the case 
of solutions (J211)-(J2S1)- In the latter case, it means that in some subspace (FigOl^a)) 
two different coexisting states have identical phase curves. Changing the external 
parameters of the system, for example the external pump frequency, we can control 
the phase structure of the coexisting states. Consequently we can make the system 
jump from a coexisting state to a periodic, quasiperiodic (beats) or chaotic state. 

6 Appendix 

The method presented is a version of the Lindstedt's method |17,(p.224), applied to 
autonomous systems in complex variables. If this method is applied to differential 
eqations of nonlinear optics written in the rotating wave approximation it gives the 
so-called c/osed solution. Otherwise, we get series solutions. Below we show how this 
method works for the system(jni)-(@I). Moreover, we present other periodic solutions 
of the selected equations of nonlinear optics. 
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6.1 Second- harmonic problem 

We find a periodic solution of Eqs.(jni)-0- Tlie problem is considered for arbitrary 
initial conditions: a(0) = a and 6(0) = (5. It is obvious that if e = the systemQ- 
(jlj) has a periodic solution given by the functions : a(t) = aexp (— itut) and 6(t) = 
/3exp (— z2u;t). Now, we suppose that if e 7^ the set of equations ©^O also has 
periodic solution with unknown frequencies f2 and 2^2. In order to avoid dealing 
with the unknown frequencies in the system ©^(Hl) we put Qt = lot . This leads to 

a(t) = « (^) = Kr) , m = ^ (^) = ^(^) ' (42) 



and 

da dAVl db dBVt 

dt dr u ^ dt dr u 



(43) 



On inserting (j42|) and (j43p into (J2I)^© we obtain the equations of motion in the 
new variables: 

— = -tujA + eA*B , 

u! dr 





^"^^ = ^2uB 0.5eA\ 
uj dr 


(44) 


In this case we can look for series solutions of the form: 






A = Ao + eAi + e'^A2 + ... 


(45) 




B = Bo + eBi + e^B2 + ... 


(46) 




n = to + ecui + e^uj2 + ... 


(47) 


On substituting (j45j) , 


PU|) and (P7() into(|44() we eet a recursive set of equations: 


io 


= -iuAo , 


(48) 


Bo 


= -i2ujBo^ 


(49) 


ii 


= -luAi Ao + A^Bo, 

UJ 


(50) 


Bi 


= -2iu}Bi-—Bo-0.5Al, 


(51) 


M 


= -iujA2 -Ai - —Ao + AlBo + A^Bi , 

CO 00 


(52) 


B2 


= -i2ujB2 - —Bi - —Bo - AoAi , 


(53) 



with the initial conditions: 



Ao{0) = a, BoiO)=P, 
Ai{Q) = , Bi{0) =0 for i > . (54) 
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The dot denotes that differentiations are with respect to r. The zero-order solutions 



are: 

Ao{t) = ae-''^^ (55) 

Bo{t) = (3e-''^\ (56) 

On substituting (|^ and (j^ into (fSUj) and (|^ we have: 

ii = -iooAi+iuJiae-''^^ + a*f3e~''^\ (57) 

Bi = -i2ujBi+t2uJi(3e-^''^^ -0.5a^e-'^'^\ (58) 

The secular terms are: iujiae~"^'^ + a*Pe~"^'^ and i2ci;i/?e~*^'^'^ — 0.5Q;^e~^*'^'^. To 
eliminate the secular terms we put 

iooia + a* (3 = , (59) 

i2uJiP-0.5a'^ = 0. (60) 

These assumptions reduce the set of equations (j37j) - (fSHj) to the form 

ii = -icjAi, (61) 

5i = -i2u;5i. (62) 

The above equations with zero initial conditions Ai(0) = and -Bi(O) = have 
trivial solutions, therefore : Ai{t) = and -Bi(t) = 0. Now, we calculate the new 
frequency il = u + euoi. From (jSUj) and (lUUj) we obtain uoi = i^ and uoi = — f ^• 

This is only possible if /5 = ±|-i/^. Therefore, we have Vt = uj±0.5e\/a*a. Finally, 
from (HSJ-d^ and ^ we get 

a{t) = ae-'^'^^^'^^^'^' + e * (63) 

b{t) = ^lJ^e-'2("^5^^'°^)* + e*0 (64) 

Now, we can consider second-order corrections. Because Ai = and Bi = 0, Eqs 
(IS21)-(ISS1) have the form: 

A2 = -iujA2 Aq = ~iuA2 + iuj2ae~"^'' , (65) 

UJ 

B2 = -i2iuB2 - —Bo = -i2ujB2 + i2u)2e~''^'^^ ■ (66) 

a; 

The secular terms are: (c<j2/co')Ao and {uj2/uj)Bq. Therefore, we have to put c<j2 = 
which leads to the following equations: 

A2 = -iujA2, (67) 

B2 = -i2ujB2. (68) 

21 



The above equations also have (with zero- initial conditions) zero-solutions: ^2(t) = 
and B2{t) = 0. Consequently: 

a{t) = ae-*^^^^^^'^)* + e * + e^ * (69) 

bit) = TW4e-*'^"^^^^'^)* + 6*0 + 6^*0 (70) 

Generally, the mathematical induction method leads to uJi = and Ai{r) = Bi^r) = 
for i > 1. Therefore, the above solutions are closed. Finally, we have: 

a{t) = ae-'("^^5v^^)* , (71) 

b(t) = ^lj^e-^2(-±l^v^^^)*. (72) 

2 V a* 

Remark.The first integral of the system 0^© is of the form ua* {t)a(t)+2ujb* {t)b{t) = 
const. The solutions of Eqs. (ffT| - (ff^ naturally implies that a*(t)a{t) = const and 
b*(t)b(t) = const (see also [T!)]). 

6.2 Another selected solutions 

The method presented also allow us to find a periodic solution if the number of 
equations is lager than two, for example |24j : 

da . ,^ 

— - = -luJaa + ecb , 
at 

db 

-— = —iubb + eca , 

at 

-— = —ioJcC — eab , (73) 

at 

where Uc = u!a + ^b- Eqs. (f7H|) are used to describe the so-called parametric optical 
processes ^E\. If the initial conditions denoted by a(0) = a, 6(0) = b and c(0) = c 
satisfy the relation c = ,^7^°f^^^ then the periodic solution is given by 

a{t) = ae V V^^^+^ 

b{t) = be y Va*a + b*bJ ^ 

C(t) = ^^a^^-i(..±ev/^^^i+F6)t_ ^^^^ 

\/a*a + b*b 

The method presented is also useful if the nonlinearities in a dynamical system 
are of different rank (for example, Kerr effect in the presence of second-harmonic 
generation ^3]): 

— = —iua + ta*b — iK,a*a? , (75) 

db , 1 r, , , 

— = -i2ujb - -ea^ . 76 
dt 2 ^ ' 
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We obtain: 

a{t) = ae-'^\ (77) 



bit) = h£«^T^yja^ + ^ le---, (7^ 



H — (Ka*a ± \/K?a*^a^ + e^^y*a 



fl = LJ . 

2 

where k = eA (A is a numerical coefficient). The assumption k = eA is necessary in 
the method presented. For A = we get the solutions (ffT|l - (ff^ . 
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